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Abstract 

The coherent stochastic resonance is observed and studied with multi-step 
periodic signal in continuous medium having two absorbing boundaries. The 
general features of this process are exihibited. The universal features at the 
resonance point are demonstrated. The kinetic behaviors around the reso- 
nance point are also presented. 
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I. INTRODUCTION 



There has been a large deal of interest in the understanding of mechanism of interplay 
between random noise and a deterministic periodic signal after the pioneering achievement 
of separation of large DNA molecules in gel medium by the application of uniform and 
time-dependent periodic electric field |]||. It has been found that with this technique large 
molecules in the size range 2 to 400 kb exihibit size-dependent mobilities. Similar ideas have 
also arisen in other types of chromatographic processes ||. 

The first passage time is a useful tool to investigate the diffusive transport property 
in a medium. The theory of first passage time has been worked out in great detail for 
both infinite medium and explicitly time-independent diffusive processes [§-§|. However, for 
explicitly time-dependent processes and in finite medium an analytic closed form expressions 
are not available. In this respect also this problem attracts much attentions to the scientific 
communities. 

The first analysis of this phenomena has been done for a random walk on a lattice 
numerically, and for a diffusive process in a continuous medium with periodic signal of small 
amplitude perturbatively 0. Their results indicate that the oscillating field can create a 
form of coherent motion capable of reducing the first passage time by a significant amount. 
This fact clearly implies that the mobility of a particle in a diffusive medium can be increased 
by the application of proper oscillating field. This phenomena thereafter is known in the 
literatures as coherent stochastic resonance(CSR). 

In order to investigate the reason for this cooperative behavior of random noise and 
deterministic periodic signal this problem has been formulated in much simpler terms by 
approximating the sinusoidal periodic signal by the telegraph signal |J and concluded incor- 
rectly that the system exihibits CSR. Subsequently it has been shown || that the telegraph 
signal can not produce CSR. It is then argued |J that the low frequency behavior could 
cause such cooperative behavior. 

In this paper we approximate the sinusoidal signal by a multi-step periodic sig- 
nal(explained below) and obtain an expression for the mean first passage time(MFPT). 
After giving the derivation of MFPT in Sec. II, the results of the calculations are discussed 
in Sec. III. First we present the general characteristics of CSR. The calculation clearly ex- 
hibits how resonance appears in our multi-step approximation and fails to show in single-step 
telegraph approximation of the periodic signal explaining the conjecture of Porra ||. The 
general characteristics of the moments in our calculation are also in agreement with the 
numerical simulation of random walk model on a lattice 0. The characteristic features of 
first passage time density function (FPTDF) for this phenomena are also presented in this 
subsection. In the next subsection we focus on the resonance point and demonstrate some 
universal features associated with it. Subsequent subsection deals with the characteristic 
changes of the physical variables as we cross, in particular, around the resonance point. 
This leads to a better understanding of this cooperative behavior. Finally, few concluding 
remarks have been added in Sec. IV. 



2 



II. DERIVATION OF THE MEAN FIRST PASSAGE TIME 



We consider diffusion in one dimension perturbed by a periodic force. The motion of the 
particle is given by the Langevin equation 

X = AsinQt + £,(t), (1) 

where X refers to the stochastic variable, A and Q are the amplitude and frequency of the 
sinusoidal signal and £(£) is a zero mean Gaussian white noise of strength D with auto- 
correlation function given by 

<i{mt')>=ZD5{t-t'). (2) 

The motion is confined between two absorbing boundaries at x = and x = L. The Fokker 
Planck equation corresponding to Eq.(P is 

^1 = -AsinQt^^- + D^fK (3) 
at ox ox 2 

with absorbing boundary conditions at x = and x = L;i.e.,p(0, t) = p(L,t) = 0. We now 
introduce the dimensionless variables 

e = (A/D)x, 9 = (A 2 /D)t, u = tt/(A 2 /D), (4) 

to write Eq.(|3]) in terms of new variables: 

—bT- = -^—dT + -505— ( 5 ) 

The boundary conditions are rewritten as p(0, 9) = p(A,6) = 0, where A = (A/D)L. In 
the following we calculate all the physical quantities in terms of these new variables and 
if required, one may translate all the interpretations in terms of the usual variables by the 
transformation equations Eq.|j. 

We next approximate the sinusoidal signal by multi-step periodic signal. The construc- 
tion is as follows. We divide the half cycle of the signal by (2p + 1) intervals so that 
each interval in the horizontal #-axis is of size (A8/(2p + 1)) with uj A 9 = it. We define 
(2p + l)numbers Sk along the vertical £ -clXIS clS 

s k = 1 2p+1 2 ^ ; k = 1, 2, ...,p (6a) 

Sp+i = 1 (6b) 
Sp+i+r = Sp+i-r ;r = 1,2, ...,p. (6c) 

Each number Sk is associated with the interval g^p- < 9 < ■£0^ with k — 1,2, (2p + 1). 
The Eq.(^|) clearly shows that 

< Si < s 2 < ... < s p < s p+ i = 1 > s p+2 > s p+2 > ... > s 2p +i > . (7) 
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Eq.([?D states that in order to reach the maximum value (= 1) of the signal from the zero 
level we have to have (p + 1) step up and from the maximum to the zero level we have 
(p + 1) step down. This is for the positive half-cycle. For the negative half-cycle similar 
constructions have been done with the replacement — > —Sk,Vk and each number — is 
associated with the interval A0[1 + ^\ < 9 < A0[1 + with k = 1,2, (2p + 1). 
This aproximation for the full one cycle of the sinusoidal signal (as shown in Figgis then 
repeated for the next successive cycles. The construction clearly shows that we get back the 
usual telegraph signal with p = 0. 

One may however note that the u> which we have defined for this approximated signal 
is not the same as that of sinusoidal signal because the Fourier transform of sinusoidal 
signal would give only one frequency while this approximated signal in the Fourier space 
corresponds to many sinusoidal frequencies specially because of its sharp discontinuities. Yet 
we urge this approximation because in each interval the equation become time- independent. 

In the future development we associate the index n for the positive half-cycle and index 
m for the negative. Index i will refer the cycle number. Since the Fokker-Planck equation 
(Eq.(|5|)) in each interval will be that for a constant bias, we can express the conditional prob- 
ability density function 9 | 9') in terms of complete orthonormal set of eigenfunctions 
w n (£) satisfying the boundary conditions u n (0) = u n (A) = 0. 

p{U\ e,9')=J2<(Ou-(e)exp[-\ n (e-9 , )\ , (8) 

n 

where 

= (2/A)*eap(±s£/2)sm-^ , (9a) 

*» = ^ + 7- < 9b > 

with s as the corresponding value of Sk in the appropriate interval where the conditional 
probability is being decomposed. The conditional probability density function in any inter- 
val, say I, can then be calculated from the previous history by convoluting it in each previous 
intervals: 

p{iu 0i I 6, 0i) = / .» / g?6-i^-2-^2 n p(&, Bj | o,-i) ■ (io) 

J J j=2 

For the negative half-cycle the calculation of probability density function is similar except 
that we have to replace the index n by m and the probability density function is decomposed 
as 

p{U\^B') = Y,^)uUOex P [-X m {B-B')\ , (11) 

m 

where the expressions for «„(£) an d A m are same as in Eqs.([|). 

The survival probability at time 9 when the particle is known to start from £ = £ at 
9 = is defined as 

s(e\£ )= /Vp(£,0Uo,o). (12) 

J 
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The first passage time density function (FPTDF) g{9) is defined as 

9(0 | 6) = ^ • (13) 

Physically, g(9)d9 gives the probability that the particle arrives at any one of the boundaries 
in the time interval 9 and 9 + d6. From this density function one can calculate various 
moments: 

< 9 j >= / d99 j g(9) . (14) 
Jo 

From Eq.(|T4]) one can easily calculate mean first passage time(MFPT) < 9 > and the 
variance a 2 =< 9 2 > — < 9 > 2 of the density function g(9). 

It is then quite straight-forward to calculate the survival probability at any interval of 
any cycle. We will write down the final formulae: 

S + (9 I f ) = C+ x exp[-X ni9 1VM (9 - 2(i - 1) A 9)} x F^iuZ ) 

+ V 1 SU/ n (2p+l)(i-l) + l ^ L n (2p+l)(!-l) + l V V / /J 1 J-V Tl(2p+l)(i_l)+l/ 

;2{i-l) A9 < 9 <[2{t-l) + —^—]A9 , (15a) 

2p + 1 

I Co) = C (2p+1) ,_ 1)+(fc+1) x e* P [-A„ (2p+1)(i _ 1)+(fc+1) (£ - 2(< - 1) A 0)] x 

ni< w« I ut >} 

L n (2p+l)(i-l) + (j + l) 1 ™(2p+l)(i-l)+j J 

i=i 

r A? 

X( 



;eX ^^2p + 1 ^^ n ( 2 P+ 1 )( i - 1 ) + ( fc + 1 ) ^ ^«(2p+l)(i-l)+(j+l))] 



Xi ? i-l('U n . (2p+1)(i _ 1)+1 ) 



[2 (*-l) + ^]A0<0<[2(;-l) + |^]A0 

fc = l,2,...,(2p-l) , (15b) 



I Co) = C (2p+1)l x ezp[-A n(2p+1)i (0 - (2. - 1) A 0)] 

xA + (w~ ,ut ) x Fi_i(«~ ) 

V re(2 p +i)i' ri( 2 p+i)(i-i)+i / * x \ n(2p+i)(i_i)+i/ 

; [2(i - 1) + A < < (2i - 1) A 9 , (15c) 



I *o) = C m(2p+1){l _ 1)+1 x exp[-A m(2p+1)(! _ 1)+1 (0 - (2. - 1) A 0)] 

x ?/ I ?y ^> 

»Tl(2p+l)(i-l) + l I "(2p+l) l 

x^4 + (w~ ,ut ) x ) 

V ™(2p+l)i' ™(2p+l)(i-l) + l^ * ■ l V rl (2p+l)(i-l) + l'' 

; (2i - 1) A 9 < 9 < \(2i — 1) H — ] A , (15d) 

2p + 1 
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S-(0 I &) = ^ (2p+1)(i _ 1)+(fc+1) x exp[-\ m{2p+m _ 1)+(k+1) (e - (2i - 1) A 61)] x 

Mm (2p+i)(i-i)+o+i) i M "i(2 P +i)(i-i)+i > ^ 

2=1 

A6? 

xexp[^ — - (fcA m(2p+1)( ._ 1)+(fc+1) - X] ^ m (2 P+ i)(i-i)+( J+ i))] 
x <^ ?/~^ I ?/~^~ 

™l(2p+l)(i-l) + l I ™(2p+l)i 
X ^ + (%2p+l)i' U n (2p+1)(i _ 1)+1 ) X -^-l( M n (2p+ i ){l _i )+ i) 

;fc = l,2,...,(2p-l) , (15e) 



S_(£ | £o) = x e*p[-A m(2p+1)i (6> - 2* A 0)] 

X >4 ( ll~ ^X <^ 7/^ I ?/^~ 

V"m (2 p + i)i' u 'm (2p+1){i _ 1)+1 ; ^ u, "i(2 P +i)( l -i)+i I «(2p+l)i 

X ^ + ( M n (2p+1)l5 M n (2p+1)(l _ 1)+1 ) X F i-l( U n (2p+1)(i _ 1)+1 ) 

■[(2i-l) + ^-^\A6 <6 <2iA6 , (15f) 



where 



C+ = / A de«+(0 , (16a) 

C m = [ A dtu m (0, (16b) 
A6? 

y4+ ( M ™ ( 2p + l) l ' M ™(2p + l)( l -l) + l) = eX P[~( 2p+ 1 )^(2 P +l)(i-l) + J 

m< W~ I M + > 

n(2 P +i)(i-i)+0+i) I n.(2p+i)(i-i)+j 

J'=l 

xexp[-(^ — -)A„ (2p+1)( ._ 1)+(j . +1) 

]} , (16c) 

A# 

^ ( M m (2p+1)l ' M m (2j , +1){l _ 1)+1 ) = e;r P[ _ (2 p+ p^"»(2 P +l)(i-l)+l] 

2p 

II{ < M m(2p+l)(i-i)+( J+ l) I M m(2p+i)(i-l) +J > 
J'=l 

A6? 

x^P[-(2^i) A -(2 P+ i) (i -i )+ o + i)]} ( 16d ) 
and the function Fj is generated through the recursion relation: 
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Fj(u„ ) = < u„ I u m > xA (uL ,u m ) 

iv n (2p+l)i+l' n (2p+l)i+l 1 m (2p+l)i v m (2p+l)i' m (2p+l)(i-l) + l > 

x < ut \ut > xA + (u~ ,ut ) 

m (2p+l)(i-l) + l 1 ™(2p+l)i v "■(2p+l)i' "(2p+l)(i-l)+l ' 

* F i-iK 3p+m - 1)+1 ) . ( 17 ) 

with F (u~J = m^(^o)- The angular bracket in any equation implies dot product of the 
corresponding functions, for e.g., 

<u+\u->= f A d£u + (Ou~(0 . (18) 





The cycle variable % runs over positive integers;i.e.,z = 1,2,3, .... The positive and negative 
symbols of the survival probabilities indicate its value over positive and negative part of the 
cycles respectively. In all these expressions, viz., Eqs.(15)- (|T7|) , any subscript either n or m 
or both whereever they appear more than once the summation over them are implied. The 
effect of history is explicit in the expressions for survival probabilities. Once the survival 
probability S{6 | £o) is obtained from these formluae, the FPTDF,MFPT and the corre- 
sponding variance are obtained by employing Eqs.(|l3|)-(|T4]). Evaluation of MFPT and other 
relevant quantities requires sum of infinite series which must be truncated in order to obtain 
a final result. Convergence of MFPT is ensured by gradually increasing the number of terms 
(i.e., number of eigenvalues) for the calculation. The process is truncated when MFPT does 
not change upto two decimal point of accuracy with the change of number of terms. 



III. RESULTS AND DISCUSSIONS 



The survival probability, mean first passage time(MFPT), corresponding variances and 
first passage time density functions(FPTDF) are calculated using the derived formulae for 
this process. The results are summarised below. 



A. General features of CSR 



The MFPT is calculated for single-step telegraph signal (p = 0) with £ = A/2. Most 
of the calculations are done with this specific value of £o- The variation of the results with 
variation of £o is also demonstrated [see the text below]. No nonmonotonous behavior is 
observed in MFPT as we vary the frequency uj. This is in complete agreement with Porra's 
observation 0. The calculation is done for the length A = 20 and the result is shown 
in the curve a of Fig.2. However, when we take p = 1, i.e., when the sinusodal signal is 
approximated by two-step periodic signal, the calculation of MFPT for the same length 
shows clearly the nonmonotonous behavior. This is shown in curve b of the same figure. 
This result clearly demonstrates that mere flipping of the bias (signal) direction periodically 
would not produce the coherent motion. As the rate of flipping increases it merely prevents 
the particle more to reach the boundaries and therefore MFPT increases monotonically. It 
may be noted that when the flipping rate is very high, the effect of signal is almost nil and 
the transport is effectively diffusive in nature. This is of course true in any type of periodic 
signal. Therefore, for any type of approximation of the sinusoidal signal or for any value 
of p, this feature would show up. In particular, for p = 1, we observe from curve b of fig.2 
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that MFPT asymptotically reaches the diffusive limit A 2 /8 (= 50 in this case). The usual 
telegraph signal offers a constant bias of maximum magnitude for the larger time than for a 
two-step approximation. Hence the particle always has a larger probability of reaching the 
boundary in short time for p = case than for p > case. Hence MFPT for p = and for 
any u is always less than for p > case. This is observed in Fig. 2. 

The application of any bias always reduces the MFPT than for the non-biased diffusion. 
In CSR we always have a competition between diffusion and oscillatory effect of the bias. 
For very large frequency as the bias effect becomes ineffective MFPT would essentially 
be guided by diffusive process. For zero frequency of the multi-step periodic signal the 
MFPT can be analytically evaluated. When it starts from the mid-point of the medium it 
expresses as < 9{uj = 0,£o = A/2) >= 0.5(A/si)iem/i(siA/4). When frequency is very small, 
the process is predominantly diffusion with constant value s\ = 0.5sin(^^) effective for 
< 9 < ^?2p+i) • However, as frequency increases slowly the probability of having increased 
bias value S2 (=1 for p = 1) before it reaches the boundary increases. This bias force reduces 
the survival probability and also MFPT. Hence one would expect a minimum to MFPT. On 
the otherhand, for usual telegraph signal (p = case), for very low frequency, from the very 
beginning bias force affects the particle with its maximum strength. When the frequency is 
very low, this constant bias diffusion continues for a longer time and there is no change-over 
of the magnitude of the bias as in the case of p = 1. After having a flip, the particle again 
suffers a constant bias diffusion in the direction opposite to the previous one. As frequency 
increases slowly, this picture remains unchanged until a stage reaches for which the flipping 
effect becomes dominant during the particle's survivality inside the medium and MFPT 
increases. This is observed in Fig. 2. 

Next we continue all our calculation with p = 2 or, with three-step telegraph signal. 
Calculation reveals that the value of MFPT does not change much from that with p = 1. 
On the otherhand, p = 2 signal approximates better than p = 1 signal. We restrict our 
calculation with p = 2 approximation of the periodic signal. 

Typical survival probability and the corresponding decay rate defined as p(9) = 
— dS }y I S(9) are plotted as function of 9 for A = 20 and cu = 0.1 in Fig.3. The plot shows 
that the survival probability [plot a] goes through plateau where the change of survival prob- 
ability is comparatively less. The decay rate p(6) [plot b] correspondingly shows a minimum 
at these points. This is a characteristic feature for CSR. This feature is in agreement with 
the numerical simulation of the process as a random walk on a lattice 0. 

Next we calculate the MFPT < 9 > and the variance o 2 as a function of frequency u> for 
different lengths(A = 10, 20, 30, 40, 50). These are presented in Fig. 4 and Fig.5 respectively. 
Both the cumulants go through a minimum as frequency rises from very low value for each 
length A. This feature is also in agreement with the lattice simulation work 0. It is 
observed that the minimum for both the moments ocuur at the same frequency for each 
length. The value of MFPT < 9 > increases with the length at all frequencies. This is 
understandable because as length increases on an average the particle will spent more time 
in the medium before reaching the boundaries. It is also observed that the frequency at 
which the minimum occurs shift towards low frequency as the length increases. It implies 
that maximum cooperation between the deterministic signal and random noise occurs at 
lower frequencies as the length increases. For low resonant frequency the particle is affected 
by the bias in a particular direction for a longer period of time before it suffers a change in 
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the direction of bias, thus more probability to cover a large distance towards the boundary 
and at this resonant frequency the probability for reaching the boundary in a short time is 
maximum because if one increases the frequency more than the resonant frequency at that 
length, the flipping rate dominates and average time taken by the particle would be more. 

Fig. 5 demonstrates the lowering of the dispersion at resonant frequencies confirming that 
the cooperation is maximum at these frequencies. Dispersion is more for higher lengths and 
as seen from the figure the dispersion merges to a specific value at very low frequency at 
various lengths. 

All the previous calculations are done when the particle starts initially from the mid 
point of the medium, i.e., £o i n Eqs.(15) is taken as A/2. At the length A = 20 the resonant 
frequency is found to be 0.1. The calculations are done one at resonant frequency and other 
two at the off- resonant frequencies (u = 0.5 and uo = 0.0) when the particle starts from 
£o = (3 A where (3 lies between and l.For zero frequency the MFPT can be analytically 
obtained. Its expression reads as < 6{uj = 0, (3 A) >= (A/si)[— (3 + (1 — exp[—Si/3A])/(l — 
exp[— S]A])]. The curves are shown in Fig. 6. It is evident that the value of < 9 > is less 
for resonant frequency (curve a) than for its value for off-resonant frequencies (curves b and 
c). As frequency increases, the maximum value of < 9 > occurs at lower values of (3 or, 
when the particle starts from the left of the interval. It is known that for pure diffusion 
the location of maximum < 9 > would occur for (3 = 0.5. Our signal starts with positive 
half-cycle and therefore the survival time of the particle would be more if the particle starts 
from the left of the interval. Of course there would be some limit, because if it starts too 
much near the left end then diffusion towards the left boundary dominates and average time 
would be less. On the otherhand, if it starts from right half of the medium, the initial surge 
of the signal helps the particle to reach the boundary more quickly. Hence average time of 
duration decreases. This fact is also in conformity with lattice simulation work [0] , although 
most of the simulations in ]7) were obtained for an uniformly distributed initial condition. 

We next calculate the FPTDF g(9) for various frequencies for A = 20 and plot the 
curves in Fig.7. The resonant frequency for this length is found to be 0.1. Before the 
resonant frequency is reached, g{9) has got two distinct peaks [Fig. 7a] and at resonance two 
peaks merge to a single large peak. After the resonance many smaller peaks in g(9) gradually 
emerge as frequency increases more than the resonant frequency [Fig. 7b]. This is a general 
characteristic of CSR. The height h and the position 9 P of the first peak as a function of uo 
are plotted in curve a of Fig.8. The figure shows that height of the first peak goes through 
a maximum as we increase the frequency while the position of that peak remains practically 
constant. The height reaches the maximum near the resonant frequency demonstrating that 
the probability of reaching the boundary is in short time is maximum near the resonant 
frequency. It is a kind of reflection of having < 9 > minimum at that frequency. Therefore 
it is a general characteristic of CSR. The height and position of the second peak before the 
resonance are drawn as curve b in Fig.8. At resonance the two peaks merge and we have 
only one peak. Just after the resonance another peak starts developing and height increases 
as frequency increases further. The position and height of the second peak after resonance 
are plotted in curve c of the same figure. The merging and the reappearance of the second 
peak is also observed as a brake or discontinuity of the dashed line in this figure. 
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B. Universal features at resonance 



In this subsection we concentrate on the behaviour of the system at the resonance point. 
We have already discussed some general characteristics of CSR in the previous subsection. 
We find that for each length, A, a corresponding frequency u* exists for which < 9 > and 
a 2 become minimum implying that the maximum cooperation between the deterministic 
periodic signal and random noise of the environment is taking place in helping the particle 
to reach the boundaries. One therefore would naturally inquire about the relation of u* 
with A. The curve of uo* as a function of A is plotted in Fig. 9. In the range of A we studied 
this curve is very well fitted with the formula 

uo* = 2/A . (19) 

The values of MFPT at resonance < 9(uo*) > is plotted against the length A in Fig. 10 
and within the range of A we consider the relation between them is fitted to 

< 9(uo*) >= 0.82A-0.14 . (20) 

Of course, there will be deviation from this linear behaviour as A decreases further because 
< 9 > can not become negative and for A = (corresponding to L = 0), < 9 > should be 
zero. 

Similarly the variance <7 2 (uo*) is plotted as a function of A in Fig. 11 and within the range 
of A we consider this curve is fitted to 

A = a[a 2 (uo*)] 2 + b , (21) 

with a = .004, b = 9.29 ± 0.82. 

We have already seen that at the resonance frequency we have one very dominant peak 
of FPTDF, (?(0)[Fig.7b]. Since it is a general feature, for each length A we should get such 
behaviour. We further observe that uo* varies inversely with A (Eq.(19)). With this fact in 
our mind when we plot g(9)/uo* as a function of [u;*(A)0], we find that curves for all lengths 
superpose over each other [Fig. 12] and the pattern of g(9)/uo* for different A or uo* is very 
similar, i.e., at particular values of [uo*9], all curves show their maxima, minima, and change 
in the behavioral patterns of the curves occur exactly at the same places of[u;*0]. Similar 
characteristics are also observed in the curves of decay rate p for different frequencies. For 
illustration we plot p as a function of [uj*9] for three different lengths(A = 20 ; curve a, 
A = 28.57; curve b, A = 15.38 ; curve c), and present in Fig. 13. Therefore it shows that 
this feature is universal and [uo*9] or the cycle number is the correct variable to describe the 
resonance behaviour. We may further note that such scaling of FPTDF would not be possible 
for any frequency other than the resonant frequencies because any frequency which is not the 
resonant frequency for one length may turn out to be the resonant frequency for some other 
length and the features of FPTDF are different for resonant and off-resonant frequencies as 
has been observed from Fig. 7a - Fig. 7b. The major dominant peaks of FPTDF g{9)/u* for 
different lengths(A = 20,28.57,35,40,44.44,50) are drawn as a function of [uj*9] in Fig.12. 
There are no overlap of these curves which can be seen as we increase the resolution. The 
lowermost curve is for A = 20 and as length increases the upper curves are generated. The 
peaks for all the curves occur nearly at a quarter of a cycle. 
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The peak height h p and full width half maximum (FWHM) for each curve is plotted as a 
function of resonant frequency. The plot is given in Fig. 14. The plot shows that except for 
very low frequency they behave linearly with the resonant frequency. 

C. Behaviour around the resonant point 

We have already seen that the cooperation between the deterministic signal and the 
random noise is maximum at the resonance point where MFPT < 9 > variance a 2 take 
minimum values and corresponding FPTDF g(9) shows a major dominant peak. What would 
happen when we change the frequency slightly above and below the resonant frequency? To 
investigate the matter we choose a particular length of the medium, A = 20. The resonance 
frequency for such length, uj* = 0.1. For this particular length we take two off-resonant 
frequencies uj = 0.07 and uj = 0.13. The curves for survival probabilities as a function 
of time 9 are plotted in Fig. 15, where the curves a,b,c are for frequencies uj = 0.1, .07, .13 
respectively. The calculation of survival probability is terminated when it takes value 1 x 10~ 3 
which corresponds to zero in our calculation. The curves clearly show that as frequency 
increases the survivality of the particle prolongs. This is quite understandable because more 
oscillations prevent the particle to reach the boundary, i.e., for higher frequency we expect 
< 9 > more. The oscillatory effect is more pronounced when time is large. For large time 
we always expect the value of S(9) more for higher frequency. This is clearly observed in 
Fig. 15. But for frequency lower than the resonant frequency MFPT < 9 > is again more. 
As MFPT is the integral of the survival probability over time, we expect a change in the 
behaviour of S(6) for lower time regime. This is shown explicitly in Fig. 16. In this figure 
we find that the survival probability is more for low frequency, uj = 0.07 (curve b) than 
for resonant frequency, uj* = 0.1 (curve a) and off-resonant frequency uj = 0.1 3 (curve c). 
Especially for curve b, the value of S is so much than that for curve a, so that area under 
the curve b is more than that for curve a. We see that near about 9 — 28, the solid curve 
crosses the dashed curve. There is only one point of crossing throughout the entire time. We 
have already argued that after this crossing point oscillatory effect of this bias dominates. 
It is then clear for low time regime diffusion process competes over the oscillatory effect. 
Again, for very low 9 the chance of having increased value of the bias in the same direction is 
more for high frequency than for low frequency. Therefore for low frequency the survivality 
is more than for the high frequency. The curves in Fig. 16 also demonstrate that. 

We have already demonstrated how < 9(uj*) > varies with A in Fig. 10. The behaviour 
is linear with respect to the length of the medium. It is of interest whether the behaviour is 
changed for off-resonant frequency. For that we choose a frequency which is not the resonant 
frequency for the length A that we consider in our calculation. The MFPT < 9{uo) > for 
that off-resonant frequency are calculated for different lengths and are plotted as a function 
of A in Fig. 17. For the range of length we consider the curve is fitted to 

A = a' < 9(uj off _ res ) > 2 +b' (22) 

with a 1 = .016 and b' = 8.68 ± .33. We may note that this particular off-resonant frequency 
would be a resonant frequency for some length, A , governed by the Eq.(19). In our case 
this off-resonant frequency corresponds to length A > 50. The curve shows that when 
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A <C Ao, the rate of change of MFPT with respect to A is more and when A approaches Ao 
or, when this frequency tends to be the resonant frequency, the rate is curbed. This could 
be a signature of approaching a coherent motion from the non-cooperative behaviour. 

IV. CONCLUDING REMARKS 

We consider a diffusive transport process perturbed by a periodic signal in continuous 
one dimensional medium having two absorbing boundaries. No perturbation of the sig- 
nal amplitude is assumed in this formulation. We showed explicitly that the cooperative 
behaviour between the deterministic periodic signal and random noise leading to coherent 
motion occurs when the time-dependent sinusoidal signal is approximated by a multistep 
periodic signal and not with single-step telegraph signal. 

Although we study the process with three-step periodic signal, the formulation is quite 
general and applicable for any approximation with arbitrary number of steps. This formu- 
lation can also be applied to any arbitrary continuous periodic signal. 

It is observed that for large time oscillation of the signal plays a dominant role in the 
transport while in the low time regime frequency dependent bias force (i.e., the chance of 
having increased values of the bias in the same direction is more for high frequency than for 
low frequency) has the key factor. For very high frequency the bias effect is practically absent 
and the motion is purely diffusive in nature. At the resonance the maximum cooperation 
between the noise and the periodic signal takes place. 

An important characteristic that we observe is that at the resonance the FPTDF for 
various lengths have similar behaviour as a function of cycle number. There is only one 
dominant peak and the peak position occurs very near to a quarter of a cycle. From Fig. 12 
we observe a slight deviation of the peak positions but we believe that if the sinusoidal signal 
is approximated by more than three-step periodic signal the position of all the peaks will be 
the same. 

There is also slight discrepancy in the position of the minimum of a 2 in comparison to 
the minima of < 6 >[Fig.4,Fig.5]. This may be due to the fact that all calculations are made 
to an end when the survival probability takes a value 1 x 10 -3 . We observe that if we cut 
off the calculations for more lower values of survival probability it does not affect MFPT 
but the variances are slightly affected. Also if one approximates the sinusoidal signal better 
than three-step periodic signal one could obtain the positions of the minima of variances at 
exactly the same places as those with MFPT. 

It is intersting to observe that the decay rate at the resonance [Fig. 13] after uj*9 = 5tt/4 
is clearly a periodic function of time. 
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FIGURES 



FIG. 1. Sinusoidal signal(dashed curve) and approximated three-step (p=2) periodic signal(solid 
curve) for the full one cycle as a function of 9. 

FIG. 2. MFPT < 9 > as a function of oj ; (a)for p=0 ;the usual telegraph signal (b) for p=l ; 
the two-step periodic signal, [A = 20, = A/2]. 

FIG. 3. (a) —lnS(9) as a function of 9 (dashed curve) (b) The decay rate, p, as a function of 9 
(solid curve), [A = 20, £ = A/2,p = 2,oj = 0.1]. 

FIG. 4. MFPT < 9{oj) > as a function of frequency oj ; (a)A = 10, (b)A = 20, (c)A = 30, 

(d) A = 40, (e)A = 50, [p = 2,£ = A/2]. 

FIG. 5. The variance a 2 as a function of oj ; (a)A = 10, (b)A = 20, (c)A = 30, (d)A = 40, 

(e) A = 50, [p = 2,£ = A/2]. 

FIG. 6. MFPT < 9 > as a function of /? for length A = 20 ; (a) for resonant frequncy oj* = 0.1 
(b) for off-resonant frequncy u; = 0.5 (c) for off-resonant frequncy oj = 0.0, [p=2]. 

FIG. 7. (a) FPTDF g{9) as a function of 9 for A = 20 before resonance for frequencies oj = 
.01, .02, .03, .07 (b) FPTDF g{9) as a function of 9 for A = 20 on and after resonance for frequencies 
oj = O.l(resonant), oj = .13,0.2,0.3 respectively, [p = 2,£o = A/2]. 

FIG. 8. Height h and the position of the peak 9 P as a function of oj ; (a) for the first peak(solid 
curve) (b) for the second peak before resonance(dashed curve) (c) for the second peak after reso- 
nance (dotted curve). 

FIG. 9. Resonant frequency oj* as a function of length A. 

FIG. 10. MFPT at the resonant frequency < 9(oj*) > as a function of length A. 

FIG. 11. The variance at the resonant frequency a 2 (oj*) as a function of the length A. 

FIG. 12. The dominant peaks of g(9)/oj* at resonant frequencies for different lengths 
(A = 20,28.57,35,40,44.44,50) as a function of oj*9. The lowermost curve is for A = 20 , and as 
length increases gradually upper curves are generated, [p = 2, £o = A/2]. 
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FIG. 13. The decay rate p for different resonant frequencies as a function of lo*9; (a) 
A = 20, uj* = 0.1 (solid curve) (b) A = 28.57, uj* = 0.07 (dashed curve) (c) A = 15.38, uj* = 0.13 
(dotted curve), [p = 2,£ = A/2]. 

FIG. 14. The height h p and full width half maximum FWHM of the peaks in FIG. 12 are plotted 
as a function of their corresponding resonant frequencies. 

FIG. 15. -lnS{6) as a function of time 9 ; (a) A = 20, uj* = 0.1 (solid curve) (b) 
A = 20, uj* = 0.07 (dashed curve) (c) A = 20, uj* = 0.13 (dotted curve), [p = 2,£ = A/2]. 

FIG. 16. S(9) as a function of 9 for the same curves as in FIG. 15. 

FIG. 17. MFPT at off-resonant frequency < 9(uj jf_ res ) > as a function of A,[p = 2,£ = A/2]. 
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